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CN ■ Abstract 



We revisited the extracting iKfc/Vcil and \V„b\ through calculating the ratios 3S{Bf -> Di-{n^,p^)/BS{B° -> Ds* D^*^^) in the perturbative 

QCD approach, which is regarded as an effective theory in dealing with non-leptonic charmed B decays. Utilizing this approach, we could 

calculate the form factors effectively as well as non-factorizable and annihilation contributions. Within the updated distribution amplitudes and 

the latest experimental data, we get | Vj,i/Vci| = 0.083 ± 0.007, which favors a bit smaller | Vj,/,| compared with the averaged PDG value, but 

,S^ • agrees well the exclusively measured values. Furthermore, we predict the branching ratio of B — > DJ p^ --^ (2.7 ± 1.2) x 10^^, which could 

I ■ be measured in B factories near future. In our calculation, the major uncertainty is due to our poor knowledge of heavy meson wave functions. 

frt*' We also comment that it is not trivial to generalize this approach to Bs system, primarily because both S'' and S, can decay into the same final 

T^^ . States D+ K^ and D+ D7 . 
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The precise measurement of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements and the three related angles, namely 
a, j3 and 7, is one of the most important topics in flavor physics uJ] . Through measurements of those observables we can test the 
standard model (SM) stringently and seek for a possible new physics beyond the SM. For instance, nonzero \Vut,\ plays critical 
role in explaining the CP violation in B and K decays. Furthermore, the ratio between \V,ii,\ and \Vch\ is very important as it is 
proportional the length of the side of the unitary triangle opposite to the well-measured angle j3 . Since the end of last century, 
the obvious example of the progress was made in measuring the ratio |V„6/Vcii| and \Vi,b\ as acutely as we can ||2[]. 

In both experimental and theoretical sides, the determinations of |V„f,| are mostly on semi-leptonic B decays, involving inclu- 
sive B — > XJv decay and exclusive B — > nlv decay. However, both of them suffer from large uncertainties due to the unknown 
non-perturbative dynamics. For inclusive B -^ XJv decay, in principle, it is the most appreciated way to extract \Vub\- In prac- 
tice, however, since \Vcb\ ^ IVu^l, the total width of B — > XJv cannot be measured due to the large charm background from 
B — > XJv decays. In most regions of phase space, where the charm background is kinematically forbidden, the hadronic physics 
enters via unknown non-perturbative functions, so-called shape functions. Although there is only one shape function at leading 
order in AQco/mb^ which can be extracted from the photon energy spectrum in B ^ XsJ, the subleading shape functions are 
modeled in the current calculations |l3l|4|]. For the exclusive decay, the determination of |y„i| requires a theoretical calculation 
of the form factors using non-perturbative methods such as Lattice QCD ||5l|6|], QCD sum rules (|7|] or light front approach ||8|], 
all of which have systematic uncertainties that are hard to be quantified. The B -^ nlv branching ratio is now well known to 6%. 
Unquenched lattice QCD calculations of the B ^ nlv form factor for q^ > 16GeV^ are available and yield ||5l|6[] 

|y„/,| = (3.4±o.2+<J:5)xio-l (1) 

Light-cone QCD sum rules are applicable for q^ < 14GeV^ and yield similar results 171]. The theoretical uncertainties in ex- 
tracting \V„b\ from inclusive and exclusive decays are very much different. In PDG ||9[], the values obtained from inclusive and 
exclusive determinations are 

\V„b\ = (4.12±0.43)x 10"^ (inclusive) (2) 

= (3.5;jJ:^) X 10"^ (exclusive). (3) 

Above two determinations are independent each other, and the dominant uncertainties are on multiplicative factors. If we weight 
the two values by their relative errors and treat the uncertainties as Gaussian, we find 19(] 

|y„fo| = (3.93 ± 0.36) X 10"-\ (combined) (4) 

whose determination is dominated by the inclusive measurement. Given the theoretical and experimental progress, it will be 
interesting to see how the inclusive and exclusive determinations develop. 

Although the complexity of QCD and weak dynamics had prevented us from measuring CKM elements through non-leptonic 
decays, many methods have been explored in determining \Vi,b\ though analyzing decays of B to exclusive non-leptonic two 
meson final states. In Refs. jldllin . the authors had pointed out that we could extract the |V„i,|/|Vfi,| and |y„f,| by studying the 
ratio of the branching fractions 

£§{B° -^ D+n-)/S§{B^ -^ D+D-). 

For decays B" -^D'fn^ andfi" -^ DfD^, eitherof decay modes involves a single isospin component / = !(/= 1/2), so the final 
state interaction is not important. Furthermore, the B" -^ D^'^^ is only induced by tree operators; for B" -^ D^D , although 



penguin type diagrams will pollute this decay, the tree operators are still dominant. In Ref. Bill] , the penguin contribution has 
been considered in the generalized factorization scheme, and the decay modes with vector meson have been also considered. 
Within contemporary data, the authors estimated the ratios and gave the theoretical upper limit of |V„/,/Vt*|- However, in the 
generahzed factorization approach, the effective color number A^^*^, a non-physical constant, should be introduced in order to 
make up the non-factorizable diagrams' contributions. In the past ten years, many efforts have been promoted to estimate 
the non-factorizabale contributions quantitatively, such as the QCD fatorization approach (QCDF) 111211 . the perturbative QCD 
approach (PQCD) il3lll4ll and the soft collinear effective theory (SCET) 11 1511 . Due to charmed meson involved, QCDF and 
SCET cannot be used here. In this letter, we reinvestigate the way proposed in 1 1111 within the PQCD approach. In Refs. jlq - 
I19I1 . the decay modes with D or D^ have been also calculated, but their results cannot be used here directly because different 
parameters have been used in different modes. Here we will adopt the consistent parameter space and the updated distribution 
amplitudes. Moreover, through this study we can cross-check the applicabihty of the PQCD with charmed meson, which is 
another starting point of this work. 



For the decay modes here we discuss, the related Hamiltonian and the corresponding Wilson coefficients are given in Ref. |12C j . 
In the PQCD approach, after the integration over momenta of positive and negative direction, the decay amplitude for B — > M2M3 
can be conceptually written as 



s^ ^ jdx\dx2dxT,b\dbib2db2bT,dbT^r 



C{t)^B{x\,b\)^M2{x2M)^Mi{xiM)H{xi,bi,t)S,{xi)e 



-Sit) 



(5) 



where bj is the conjugate space coordinate of intrinsical transverse momentum of light quark, and x,- is momenta fraction of 
light quark in each meson. Tr denotes the trace over Dirac and color indices, C{t) is the Wilson coefficient evaluated at scale t, 
and H{xi,bi,t) is the hard QCD part, which can be calculated perturbatively. The function 4>m is the wave function of a meson 
M, and the function St{xi) describes the threshold re-summation which smears the end-point singularities on x,. The last term, 
e^^^'\ is the Sudakov form factor which suppresses the soft dynamics effectively. More detailed information on the PQCD can 
be refereed to Ref. 11411 . In the PQCD picture, the decay amplitudes of related decay modes are written as 

^(fiO ^ d;d+) = -^{y,,x[^^i^(C2 + ^Ci) +F^,^(Ci)] - y,iy,:.[7f ^(C4 + ^Cs +Cio + ^Cg) +F,^f (C3 +C9) 

+f/^(C6 + ^Cs + Cg + ^Ct) + F,f (Cs + C7) + -F,f^(C4 + ic3 - ^Cio - ^Cg) 
3 3 3/0 



-^if (C3 - \c,) + F,fiCe + ^Q - ^Q - ^Cv) + FJ;f (C5 - ^Cj)]}, 



(6) 



^(BO^D77r+) = -^y„,y;{F,^^(C2 + ici) + F^^(Ci)}. (7) 

In above formulae, F^,\ is the amplitude from (non)-factorizable emission-diagrams with {V —A){V —A) (or {V —A){V +A)) 
operators. SP stands for the operators —2{S — P){S + P), which comes from the Fierz translation of {V —A){V +A) of different 
color structure. Also, a{n) means the contribution from (non)-factorizable annihilation diagrams. The detailed expression of 



each F can be found in Refs. 117 



1911 . and we would not list them explicitly here. Note that the /v's are different for different 
decay modes, as they are mode-dependent. In the calculation, the parameters we input are the non-perturbative parts, namely the 
distribution amplitudes and decay constants. For the light mesons n and p, the distribution amplitudes are well studied within the 
light-cone QCD sum rules, and we will use the latest versions 12111 . However, for heavy mesons B and D, their wave functions 
and distribution amplitudes are still less known, even though much efforts have been made 1221 . For the B meson distribution 



TABLE I: The numerical results of T and P for each decay mode. Theoretical Results are the predictions of the branching ratios in the PQCD 
approach. For comparison, we also list the experimental data. 



Mode 


T 


P 


Theoretical Results 


Experimental Data \9} 


B"^Z)7;r+ 


1.56-0.01 i 





(2.6+!!:^+ii) X 10-^ 


(2.4±0.5)x IQ-^ 


B" -^D-p+ 


1.64-0.04/ 





(3.1+i;°+|:,^)xlO-^ 


<2.4x 10-4 


B^^D*-K+ 


1.47-0.02/ 





(2.7+!jn^)xl0-^ 


(2.6+l]i) X 10-5 


fi" -^ D-D+ 


2.67-0.17/ 


0.23-0.01 / 


(7.7+|:^l«) ^ 10-3 


(7.4±0.7)x 10-3 


fiO -^ D-D*+ 


2.68-0.17/ 


0.23 - 0.02 / 


(7.6+!::+«)xio-3 


(8.2±l.l)x 10-3 


S" -^ D*-D+ 


2.72-0.21 / 


0.25 - 0.02 / 


{lJt\i+_tj) X 10-' 


(7.6±1.6)x 10-3 



amplitude, we adopt the model 11 1311 : 

(/)B (x, b) = Nbx^ ( 1 - x)2 exp 



1 ^xMbY colb^ 



(8) 



n 



with the shape parameter cob ~ 0.40 ±0.05 GeV, which has been tested in many channels such as B ^ rnz^Kn 11411 . The 
normalization constant A^B is related to the decay constant /b = 190 Me V [13]. As forDorDy meson, the distribution amplitude, 
determined in Ref. 11711 by fitting, is 



VD = 



2 76 



/d6x(1 -x) [1 + Cd( 1 -2x)]exp 



aP-b^ 



(9) 



where Cd = 0.5,0 = 0.1 forD meson and Cd = 0.4,0) = 0.2 forDj. 
Because of the unitarity of the CKM matrix, we have a relation: 



(10) 



Within the parametrization proposed by Wolfenstein, IVHfcVi*^! is about X'^, while |VfftV,*| and |VtiK*sl ^""e about A^, where 
X ~ 0.22. Therefore, |V„/:,V,*i| can be ignored safely, which leads to VtbV*^ = —VcbV*^. For simplicity, Eqs. (|6H3 could be 
re-written as; 



£/{B' 



,0 



D;D+) = V,hV:,{TDD,+PDDj, 



(11) 

(12) 



7i)D,, PdDs and Td.ji can be calculated directly in the PQCD approach, whose central values are listed in Table. 1. In Table 



21, i.e. \Vub\ = 



1, we also list the branching ratios calculated within the PQCD approach with the parameter values 
0.00359, \Vcb\ = 0.041 and |y„| = 0.973. The first error is from the threshold re-summation, and the second one comes from 
the distribution amplitude of heavy meson. For uncertainties of heavy meson, please look at Refs. lll7l - ll9ll . Comparing our 
theoretical predictions with experimental data, we find that our predictions are quite close to the data after combining all kinds 
of uncertainties, even though the central values are a bit larger than those of experimental data, which may hint a smaller \Vub \ ■ It 
should be also noted that our result of B" -^ D^ 7r+ is smaller than that from Ref. llisll . mainly because the different light meson 
distribution amplitudes are used. 



Q, 



Following Ref. Illlll . we define the ratios and estimate the values as. 



R 



n/D 



R 



p/D 



R 



7l/D* 



R 



P/D- 



R 



n/D - 






V, 



lib 



Vet 
Vuh 



Td,% 



K 



cb 



TdsD + PdsD 
Td,p 



V, 



lib 



Vet 
Vut 



Td,d + Pd,d 



Vet 



Td,d* + Pd,d* 
Td,p 



V, 



lib 



Vet 



TdsD* + PdsD 

Td*j: 



El 

p? 



p? 



PL 
p?* 

2 / p 
Pc 

P?* 



Vub 



Vet 



yub 



Vet 



X (0.44 ±0.07), 



X (0.53 ±0.09), 



V 



lit 



Vet 

Viit 



X (0.46 ±0.08), 



El 
P? 



V 



Vet 

2 



lit 



Veb 



X (0.55 ±0.09), 
X (0.45 ±0.07), 



(13) 
(14) 
(15) 
(16) 
(17) 



Td*d + Pd*d 

where p^ is the cm. momentum of the decay particle X. Combining the experimental data of Table 1 , we arrive at the following 
ratios: 



-2. 



^;i/D = (0.32 ± 0.07) xlO"-"; 
R„/D* = (0.29 ±0.07) X IQ-^; 

^;i/D = (0-34±0.10)xlO"'. 

After compared with the predictions of the PQCD approach ( fT3]fT7] i. we can get: 



(18) 
(19) 
(20) 



V„b 



Vet 



0.085 ±0.011 from R^/d 



= < 0.079 ±0.011 from R 



k/D'- 



(21) 



0.087 ±0.014 from R^/u 



and the averaged value is 



Vub 



Vet 



0.083 ±0.007. 



If we set \Vcb\ =0.0415 ±0.0011 which is the global fit PDG result, we can get 

|V„i| = (3.44±0.30)xlO"l 



(22) 



(23) 



Compared with the PDG values Eqs. (|2]-01i, our central value favors a bit smaller \Vut\, but agrees well with the exclusive data. 
We also notice that the uncertainty is still large, which is because there are larger uncertainties in both theoretical calculations 
and experimental data. With much more knowledge in QCD and precisely measured experimental data in future, the errors in 
our prediction will be reduced. 

With the averaged ratio of |V„i/ycf,| Eq. ( |22] |. by using both a relation (fT4l i and experimental data for B" — > D^D^, we can 
predict the branching ratio of B*' -^ DJP^ as: 



^B° ^ D; p+) ^ {2.1 ±1.2) X 10 



-5 



(24) 



which would be measured easily in future B factories. Please also note that our theoretical estimate within the PQCD approach 
for the branching ratio, as shown in Table 1, is 



?(B''^D;p + ) = (3.1+',:;;+|;^)xlO 



(25) 



To our best knowledge about this channel, still there are only upper limits from CLEO 12411 and ARGUS 12511 . We hope future B 
factories could measure this channel to test our prediction. 

Finally we add one comment on By decay: With the help from CDF, D0 and LHC, the measurement of branching ratios of Bj 
decay modes becomes much more precise than before. In principle, we can generalize the strategy of B^/ decays to By by using 
SU(3) symmetry. We can define the ratio similarly: 

However, both decay modes B^ -^ D^K^ and B^ -^ DjDf are much more complicate than decay modes of B^. Primarily 
because both Bj and B^ can decay to the same final state DjK^, we cannot discriminate whether the final states come from Bj 
or By. The same situation occurs in By — > D^Df decay mode. Furthermore, the By -^ D^K^ mode also has the annihilation 
diagrams, which will lead to larger uncertainty in calculating the R^j^ij)- In practice, in order to measure the CKM angle 7, the 
experimentalists often analyze the time dependent CP violation of these kinds of decay modes. 

As summary, we revisited extracting IVuii/Vcf,! and \Vub\ through calculating the ratios ^(B" — > D^ {^^ ,P^)/ ^{B^ -^ 
Dl* d'^*') in the PQCD approach. Using the updated distribution amplitudes and the latest experimental data, we got 
\Vub/Vcb\ = 0.083 ± 0.007, which favors a bit smaller |V„fc| compared to the value of PDG. Moreover, we predicted the branch- 
ing ratio of B" -^ Djp^ ^ (2.7 ± 1.2) x 10^^, which can be measured in future B factories. Furthermore, we could test the 
applicability of the PQCD with charmed mesons through comparing the future data and our calculation. 
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